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As has been confirmed meanwhile by lattice-QCD calculations (see e.g. Ref. [1]), the 
confinement spectrum of non-exotic quark-antiquark systems has its ground state for scalar 
mesons well above 1 GeV in the Resonance Spectrum Expansion (RSE) 1 . For instance, 
in the 5-wave Kir RSE amplitude, a broad resonance was predicted slightly above 1.4 
GeV [2], which is confirmed by experiment as the Kq(1430). However, a complete nonet 
of light scalar mesons was predicted [3] as well, when a model strongly related to the RSE 
and initially developed to describe the cc and bb resonance spectra [4] was applied in the 
light-quark sector. Thus, it was found that the light scalar- meson nonet constitutes part of 
the ordinary meson spectrum, albeit represented by "extraordinary" [5] poles [2]. Similar 
resonances and bound states appear in the charmed sector [6], and are predicted in the 
S-meson spectrum [7,8]. 

A recent work [9] confirmed the presence of light scalar-meson poles in the RSE am- 
plitude for 5"-wave and P-wave irir and Kir contributions to three-body decay processes 
measured by the BES, E791 and FOCUS collaborations. 

1 Scattering poles 

It is generally accepted that resonances in scattering are represented by poles in the "second" 
Riemann sheet of the complex energy plane [10]. Let us assume here that in a process of elastic 
and non-exotic meson-meson scattering one obtains scattering poles at 



1 The RSE was designed for the description of the complete resonance structure in meson-meson scattering, for 
both the heavy- and light-quark sectors. 
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Simple poles in S may be considered simple zeros in its denominator. Hence, assuming a poly- 
nomial expansion, we may [11,12] represent the denominator D of S by 



(2) 



D(E) oc (E-P )(E-P 1 )(E-P 2 )... . 

Unitarity then requires that the S"-matrix be given by 2 

(E - P*) (E - P*) (E - P*) . . . 

S(E) = . (3) 

(E-P ) (E-P l )(E-P 2 )... 

If we assume that the resonances (1) stem from an underlying confinement spectrum, given 
by the real quantities 

E = E , E u E 2 , ... , (4) 

then we may represent the differences (P n — E n ), for n = 0, 1, 2, . . ., by AE n . Thus, we obtain 
for the unitary S'-matrix the expression 

(E-E - AE *) (E -Ex- AEx) (E — E 2 — AE 2 *) . . . 

S(E) = . (5) 

(E-E - AE Q ) (E-Ex- AEx) (E-E 2 - AE 2 ) . . . 

So we assume here that resonances occur in scattering because the two-meson system couples 
to confined states, usually of the qq type, viz. in non-exotic meson-meson scattering. Let the 
strength of the coupling be given by A. For vanishing A, we presume that the widths and real 
shifts of the resonances also vanish (see Fig. 1). Consequently, the scattering poles end up at the 
positions of the confinement spectrum (4), and so 



AE n 



for n 



0, 1, 2, 



(6) 



A|0 



As a result, the scattering matrix tends to unity, as expected in case there is no interaction. 
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Figure 1: Illustration of a scattering pole and the related energy level Eq of the confinement spec- 
trum (4). 

An obvious candidate for an expression of the form (5) looks like 

1 + A2 



E — En 



S(E) 



1 + A 2 IE 

I n 



G(E) 



(7) 



2 Note that we do not consider here a possible overall phase factor representing a background and stemming 
from the proportionality constant in formula (2). 
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where G is a smooth complex function of energy E, and where, at least for small values of the 
coupling constant A, one has 

AEi « -X 2 G (E.i) for i = 0, 1, 2, . . . . (8) 

Relation (8) can be easily understood, if we assume that for small A poles show up in the vicinity 
of the energy values (4) of the confinement spectrum. As a consequence, at the zero Pi of the 
denominator, near the i-th recurrency of the confinement spectrum Ei, the term n — % dominates 
the summations in formula (7), i.e., 




For larger values of A, one cannot perform the approximation Pj « E{ in Eq. (9). In such cases, 
the left-hand part of Eq. (9) must be solved by other methods, usually numerically. However, 
since it is reasonable to assume that poles move smoothly in the lower half of the complex energy 
plane as A varies, we may suppose that the left-hand part of Eq. (9) has solutions which, when the 
value of A 2 is continuously decreased, each correspond to one of the values out of the confinement 
spectrum (4). 

When all scattering poles in expression (5) are known, one can — with unlimited accuracy - 
determine the function G in formula (7). Once G is known, one can search for poles by solving 
the left-hand part of Eq. (9). However, further restrictions can be imposed upon expression (7). 
For a two-meson system, there may exist bound states below the meson-meson scattering thresh- 
old. Such states are represented by poles in the analytic continuation of expression (7) to below 
threshold, on the real axis in the complex energy plane. Consequently, in the case that a con- 
finement state, say Eq, comes out below threshold, its corresponding pole is, at least for small 
coupling, expected to be found on the real axis in the complex energy plane. Using formula (8), 
we obtain 

G {Eq) real for E < threshold. (10) 

Moreover, in order to ensure that scattering poles come out in the lower-half of the complex 
energy plane, also using formula (8), we find that above threshold G must be complex, with a 
positive imaginary part. 

2 Partial waves 

In different partial waves, resonances come out at different masses. At threshold, where the 
total invariant mass of the two-meson system equals the sum of the two meson masses, one has 
additional conditions. For S waves, since cross sections are finite, we must demand that G do 
not vanish at threshold, whereas, for P and higher waves, as cross sections do vanish, G should 
vanish as well. 

A possible expression that satisfies all imposed conditions reads 

ipaje{pa) h^\pa) for £ = 0, 1, 2, . . . , (11) 

where p represents the linear momentum in the two-meson system and a a scale parameter with 
the dimensions of a distance. The well-known scattering solutions ji and h!p stand for the 
spherical Bessel function and the Hankel function of the first kind, respectively. 
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Thus, we arrive at a good candidate for a scattering amplitude of resonant scattering off a 
confinement spectrum, reading 





Upajjipa) 



(12) 



l + 2zA 2 1^^^-) 
I n E — E n ) 



Hpaj e (pa) h^ipa) 



where we have introduced the two-meson reduced mass /i and, moreover, relative couplings g n g, 
which may be different for different recurrencies of the confinement spectrum. 

As it is written, formula (12) seems to allow a lot of freedom, through adjustments of the g n (_ 
to experiment. In principle, it might even be useful to carry out such data fitting, so as to gain 
more insight into the details of the coupling between a two-meson system and a confined qq state. 
However, experimental results are so far much too incomplete to make a detailed comparison to 
our expression possible. 

The spin structure of quarks, besides being important for the spectrum of a qq system, is 
also crucial for the short-distance dynamics, hence for the properties of the coupling between 
qq and meson-meson states. In the 3 Pq model [13,14], it is assumed that a two-meson system 
couples to a qq state via the creation or annihilation of a new qq pair, with vacuum quantum 
numbers J PC = ++ . Under this assumption, all relative couplings can be determined from 
convolution integrals of the wave functions. In Refs. [15,16], such integrals have been calculated 
for general quantum numbers, including flavour. The latter results leave no freedom for the 
coupling constants in formula (12), except for an overall strength A, which parametrises the 
probability of qq creation/annihilation. 

This way, the full spin structure of the two-meson system is entirely contained in the relative 
coupling constants g n £. Yet, direct comparison of the results given in Refs. [15,16] to experiment 
would still be of great interest. 

The relevant qq <-> MM coupling-constant book-keeping has been developed in Refs. [15,16]. 
The latter scheme not only eliminates any freedom, but also — by construction — restricts the 
number of possible MM channels that couple to a given qq system. Nonetheless, the number of 
involved channels rapidly grows for higher radial and angular excitations of the qq system. 



The scattering matrix is not directly observable, but only through quantities like cross sections 
and production rates. It is straightforward to determine cross sections [17] and, after some 
algebra, production rates [18] from expression (12). However, a complete modelling of strong 
interactions is more complex. For example, a cc vector state couples, via OZI-allowed decay, to 
DD, but also to DD*, D*D*, D S D S , . . . [4]. Consequently, the involved two-meson channels 
couple to one another as well. So the first extension necessary for a more proper description of 
strong interactions is the formulation of a multichannel equivalent of expression (12). This issue 
has been dealt with in Ref. [19]. It involves coupling constants similar to the ones discussed 
above, but now for each two-meson channel. 

A meson-meson channel is characterised by quantum numbers, including flavour and isospin, 
and the meson masses. However, many of the needed masses are unknown yet, while most mesons 
only exist as resonances. 
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Observables 
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In experiment, one can concentrate on one specific channel. On the other hand, in a mean- 
ingful analysis all channels that couple must be taken into account. For example, one may argue 
that for the description of tttt scattering below the KK threshold the channels KK, 7777, . . . can 
be neglected. But then one ignores virtual two-meson channels, which may have a noticeable 
influence below the KK threshold. 

Furthermore, qq states may couple to one another via meson loops. Typical examples are: cc 
vector states, which become mixtures of 3 Si and 3 D\ via loops of charmed mesons, and isoscalar 
qq states, where kaon loops mix the uu + dd and ss components. One then obtains different 
interplaying confinement spectra, which may become visible in production rates. The extension of 
expression (7) to more than one qq channel has been considered in Refs. [3,20], for the description 
of the a and /o(980) resonances. 

4 The parameters 

Besides the parameters A and a, formula (12) contains an infinite number of parameters E n %. 
These represent the unknown and even hypothetical spectra of confined qq systems. From exper- 
iment, we only have data at our disposal for resonances in meson-meson scattering or production. 
Formulae like expression (12) are intended to interpolate between the observed resonances and 
the underlying — largely unknown — confinement spectrum. 

In Fig. 2 of Ref. [21] (see Fig. 3), we showed, for 5-wave isodoublet Kn scattering, how cross 
sections determined by the use of formula (12) vary with increasing values of the coupling A. For 
small A, the nonstrange-strange (ns) confinement spectrum is well visible in the latter figure, 
whereas for the model value of the qq <->• meson-meson coupling experiment is reproduced. 

A 2 ^ 

.33 .67 1.0 exp 

1=1 -0(4415) 

1 1 ^(4160) 

IZZZZZZ=ZZ=Z iz=] ^(4040) 

— ~ZZT~ZZT~.ZZ~ i ^(3770) 

^(3685) 



— — w 

J PC = 1 , cc states 

Figure 2: The theoretical values of the central resonance positions for charmonium S and D states for 
various values of the model parameter A, compared to the experimental situation. The various dashed 
lines indicate the threshold positions of the strong decay channels DD, DD*, D S D S , D*D*, D S D* and 
D*D*. 

Furthermore, in Fig. 3 of Ref. [22] (see Fig. 2) we showed a similar behaviour as a function 
of A for J PC = 1 cc states. For A = 0, we find the theoretical ground state at 3.46 GeV, 
whereas for A = 1 it coincides with the experimentally observed J / ip mass. The model employed 
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to determine the results of this figure was a multichannel extension of formula (12), taking 
moreover into account the degeneracy of certain confined qq states. 

From these results we may conclude that, although there is some connection between the 
confinement spectrum (A = 0) and the resonances and bound states of two-meson systems (A = 
1), it is not a simple one-to-one relation. Moreover, the level splittings of the confinement 
spectrum appear distorted in experiment. In particular, the experimental ground states show up 
much below the ground states of the hypothetical confinement spectrum. 

Over the past decades, many models have been developed for the description of meson spectra. 
Only very few of those models are based on expressions for two-meson scattering or production. 
Here, it is stressed that no data for the spectra of confined qq systems exist. We only dispose 
of data for resonances in meson-meson scattering or production [3,4,17]. Nevertheless, in order 
to unravel the characteristics of the qq confinement spectrum, we must rely on results from 
experiment, even though the available data [23] are manifestly insufficient as hard evidence. 

We observe from data that the average level splitting in cc and bb systems equals 350-400 
MeV, when the ground states, J ftp, rj c , T(IS') and rjb are not taken into account [4]. Furthermore, 
mass differences in the positive-parity /2 meson spectrum, which are shown in Table 3 of Ref . [24] , 
hint at level splittings of a similar size in the light qq spectrum. In Ref. [24], possible internal 
flavour and orbital quantum numbers for /2 states were discussed. 

Moreover, the few available mass differences for higher recurrencies indicate that level split- 
tings might turn out to be almost constant for states higher up in the qq spectra as well [23,25], 
a property shared by the spectrum of a simple non-relativistic harmonic oscillator. Over the 
past thirty years, we have systematically discussed an ansatz for harmonic-oscillator confine- 
ment. A formalism which naturally leads to a harmonic-oscillator-like qq confinement spectrum 
starting from QCD, by exploiting the latter theory's Weyl-conformal symmetry, can be found in 
Refs. [26,27]. 

Guided by the — not overwhelmingly compelling — empirical evidence that level splittings 
may be constant and independent of flavour, and given the obvious need to further reduce the 
parameter freedom in expression (12) for the two-meson elastic T-matrix, we simply choose here 
the qq level splittings E( n+ i)i — E n i to be constant and equal to 380 MeV, for all possible qq 
flavour combinations. The remaining set of parameters Eqo, different for each possible qq flavour 
combination, can be further reduced [17], via the choice of effective valence flavour masses and 
a universal frequency uj. In the future, when more data become available on the spectra of 
qq systems, higher-order corrections to the harmonic-oscillator spectrum may be inferred. At 
present, this does not seem to be feasible. 
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5 S'-wave scattering for J = 1/2 
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Figure 3: Cross section for S- wave isodoublet Kit scattering. Left: For very small values of A, one 
observes the J PC = ++ ns confinement spectrum. Middle: When A takes about half its model value, 
one notices some more structure for low invariant masses. Right: At the model's value of A, this 
structure is dominant and well in agreement with the experimental observations. The data are taken 
from Ref. [28] (open circles) and Ref. [29] (full circles). 

In Fig. 2 of Ref. [21] (see Fig. 3), we compared the result of formula (12) to the data of 
Refs. [28,29]. We observed a fair agreement for total invariant masses up to 1.6 GeV. However, 
one should bear in mind that the LASS data must have larger error bars for energies above 
1.5 GeV than suggested in Ref. [29], since most data points fall well outside the Argand circle. 
Hence, for higher energies, the model should better not follow the data too precisely. 
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Figure 4: S- wave Kr] "cross section" (see text), as a function of total invariant mass, (a): From 
threshold up to 2.1 GeV. (b): Detail for lower energy, (c): Detail for higher energy. 



Now, in order to have some idea about the performance of formula (12) for 1=1/2 S-wave 
irK scattering, we argue that, since in our model there is only one non-trivial eigenphase shift 
for the coupled 7rK+r)K+r] ; K system, we may compare the phase shifts of our model for t]K 
and rj'K to the experimental phase shifts for Kn. We did this comparison in Figs. 6 and 7 of 
Ref. [21] (see respectively Figs. 4 and 5), where, instead of the phase shifts, we plotted the cross 
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sections, assuming no inelasticity in either case. The latter assumption is, of course, a long shot. 
Nevertheless, we observe an extremely good agreement. 
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Figure 5: S-wave Kr/ "cross section" (see text), as a function of total invariant mass. 

Apparently, we may conclude that the phase motion in the coupled irK, t]K and rj'K system 
is well reproduced by the model. In particular, one could have anticipated that irK and r]K have 
very similar phase motions, because rjK has been observed to almost decouple from irK. This 
implies that the corresponding T-matrix for a coupled ttK+i^K system is practically diagonal. 
Knowing, moreover, that this system has only one non-trivial eigenphase, we should then also 
find almost the same phase motion for irK and r]K. 

6 The 7 = 1/2 S- wave poles 

Since the model reproduces fairly well the data for the 1 = 1/2 S'-wave, it is justified to study 
its poles. In Table 1 we collect the five lowest zeros of formula (9). Only three of the five 



Pole Position (GeV) 


Origin 


0.77- 0.28i 


continuum 


1.52 - O.lOi 


confinement 


1.79 - 0.05i 


confinement 


2.04-0.15z 


continuum 


2.14-0.07i 


confinement 



Table 1: The five lowest zeros of formula (9). 

corresponding poles are anticipated from the J p = + ns confinement spectrum, coming out at 
1.39 GeV, 1.77 GeV, 2.15 GeV, ... . So we expected only three, but find five poles in the invariant- 
mass region below 2.2 GeV. This shows that the transition from formula (5) to formula (12), is 
not completely trivial. A forteriori, expression (9) even has more zeros than expression (2). It is 
amusing that Nature seems to agree with the form of the scattering matrix in formula (12). As a 
matter of fact, the latter expression can be obtained by a model for confinement [4,30], whereas 
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formula (5) only expresses one of the many possible ways to obtain poles in the scattering matrix 
at the positions (1). 

The extra poles (continuum poles), which disappear towards negative imaginary infinity when 
the overall coupling A is switched off, can be observed in the experimental signal by noticing the 
shoulders at about 1.4 GeV in ttK scattering (see Fig. 2 of Ref. [21] (Fig. 3)), and at about 1.9 
GeV in rj'K (see Figs. 6 and 7 of Ref. [21] (respectively Figs. 4 and 5)). The shoulder in irK 
corresponds to the confinement state at 1.39 GeV, on top of the larger and broader bump of the 
continuum pole at (0.77 — 0.28z) GeV, while the shoulder in rj'K corresponds to the continuum 
pole at (2.04 — 0.15i) GeV, on top of the larger and broader bump of the confinement state at 1.77 
GeV. Such subtleties in the data may have been overlooked in the corresponding Breit-Wigner 
analyses. 

There is one more observation to be made at this stage. The central resonance peak of the 
lower enhancement in S"-wave irK scattering (see Fig. 2 of Ref. [21] (Fig. 3)) is at about 830 
MeV, whereas the real part of the associated pole is at 772 MeV. Hence, identifying the real part 
of the pole position with the central peak of a resonance may be quite inaccurate. 

With respect to the positions of the poles given in Table 1, it must be stressed again that 
these are model dependent. So the model (12) only indicates the existence of such poles in the 
respective regions of total invariant mass. A more sophisticated model, which fits the data even 
better, will find the poles at somewhat different positions. 
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